1/1 


'HD-A127  776 
UNCLASSIFIED 


RESEARCH  IN  ALGEBRAIC  MANIPULATIONS)  MASSACHUSETTS 
INST  OF  TECH  CAMBRIDGE  LAB  FOR  COMPUTER  SCIENCE 
J  MOSES  14  DEC  S3  AFOSR-TR-84-0027  AFOSR-80-0250 

F/G  12/1  NL 


AFOSR-TR. 


security  classification  of  this  page 


i«  report  security  classification 

>  UNCLASSIFIED 


2a  security  classification  authority 


2b.  oeclassification/oowngraoing  schedule 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBERtS) 


:  7 


REPORT  DOCUMENTATION  PAGE 


tb.  restrictive  markings 


b.  OFFICE  SYMBOL 
(If  applicable) 


NAME  OF  PER 

assachusetts  Institute  of 
echnoloev 


6c.  AOORESS  (City.  State  and  ZIP  Coda) 

Laboratory  for  Computer  Science,  545  Tech¬ 
nology  Square,-  Cambridge  MA  02139 


8b.  OFFICE  SYMBOL 
(If  applicable) 


12.  PERSONAL  authorisi 

Joel  Moses 


13a  type  of  REPORT 

Interim 


16.  supplementary  notation 


3.  distribution/availability  of  report 

Approved  for  public  release;  distribution 
unlimited. 


5.  MONITORING  ORGANIZATION  REPORT  NUMBERISI 

AFOSR-TR.  7 


?«.  NAME  OF  MONITORING  ORGANIZATION 

Air  Force  Office  of  Scientific  Research 


?b.  ADDRESS  (City.  Slate  and  ZIP  Code) 

Directorate  of  Mathematical  &  Information 
Sciences,  Bolling  AFB  DC  20332 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

AF0SR-80-0250 


10.  SOURCE  OF  FUNOING  NOS. 


PROGRAM 

PROJECT 

ELEMENT  NO. 

NO. 

61102F 

2304 

1? 

COSAT  I  COOES  i 

PIE  LD 

GROUP 

SUB  GR 

— 

-  - 

13b.  TIME  COVERED  14.  DATE  OF  REPORT  (Yr.  Mo..  Day) 

FROM  1/7/82  TO 30/6/83  14  DEC  83 


18.  SUBJECT  TERMS  ^Continue  on  reverse  if  neceeeary  and  identify  byt 


ms 


19  ABSTRACT  'Continue  on  reverse  if  necessary  and  identify  by  block  number/ 

Professor  Shunro  Watanabe  of  'span  has  been  visiting  the  group  since  March  1982.  In  Ithe 
past  year  he  has  been  able,  with  the  support  of  the  group,  to  develop  a  major  new  program 
for  solving  in  closed  form  a  wide  variety  of  ordinary  differential  equations. 

The  MACSYMA  system  uses  two  large  packages  for  solving  ODEs.  The  first  package  was  written 
by  Jeffrey  Golden  of  the  group.  It  solves  equations  by  using  a  few  well  known  techniques 
(e.g.,  separation  of  variables).  It  is  quite  weak  in  dealing  with  second  order  equations, 
such  as  hypergeometric  functions.  The  second  package,  due  to  Ed  Lafferty  of  the  MITRE 
Corporation,  incorporates  many  more  techniques,  but  is  still  not  a  very  general  solver  of 
ODEs. 


Professor  Watanabe' s  approach  is  to  convert  most  second  order  ODEs  into  (CONTINUED) 


20  DISTRIBUTION,  AVAILABILITY  OF  ABSTRACT  21  ABSTRACT  SEC 

UNCLASSIFIEO/UNLIMITEO  E  SAME  AS  RPT  KdtiCUSERsD  MI  1C  L/\.J  1  I  ED 

22*.  NAME  OF  RESPONSIBLE  INDIVIDUAL  22b  TELEPHONE  NUMBER 

t include  .4nra  C  ode  i 

Dr.  Robert  N.  Buchal  <  ■(>.  707-  4939 


21  ABSTRACT  SECURITY  CLASSIFICATION 


DD  FORM  1473,  83  APR 


EDITION  OF  1  JAN  73  IS  OBSOLETE 


|22l  office  symbol 


- 

I'M  '.LASSIMCATION  Of  THIS  PAG1 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  FACE 


ITEM  #19,  ABSTRACT,  CONTINUED:  instances  of  the  so-called  P-function.  These  functions, 
originally  due  to  Riemann,  have  been  extensively  studied  by  the  Japanese  mathematician 
Fukuhara  (al-so  spelled  Hukuhara).  For  the  past  year,  Watanabe  has  implemented  a  package 
that  transforms  a  large  class  of  equations  into  P-functions.  He  has  used  Kame's  table 

as  a  test  bed. 


Accession  For 

NTIS  GRAAI 
DTIC  TAB 
Unannounced 
Justification. 


By - — — - 

Distribution/ 

Availability  Codes 
Avail  and/or 
Dist  Special 


UNn.ASSIUlT' 


AFOSR-TR- 


;  A 


^  -j  7 

PROGRESS  REPORT  1982/83 

Professor  Shunro  Watanabe  of  Japan  has  been  visiting  our  group  since  March, 
1982.  In  the  past  year  he  has  been  able,  with  the  support  of  our  group,  to 
develop  a  major  new  program  for  solving  in  closed  form  a  wide  variety  of 
ordinary  differential  equations. 

-Si 

The  MACSYMA  system  uses  two  large  packages  for  solving  ODEs,  The  first 
package  was  written  by  Jeffrey  Golden  of  our  group.  It  solves  equations 
by  using  a  few  well  known  techniques  (e.g.,  separation  of  variables).  It 
is  quite  weak  in  dealing  with  second  order  equations,  such  as  hypergeometric 
functions.  The  second  package,  due  to  Ed  Lafferty  of  the  MITRE  Corporation, 
incorporates  many  more  techniques,  but  is  still  not  a  very  general  solver  of 
ODEs. 

V' 

Professor  Watanabe's  approach  is  to  convert  most  second  order  ODEs  into 
instances  of  the  so-called  P-function.  These  functions,  originally  due  to 
Riemann,  have  been  extensively  studied  by  the  Japanese  mathematician 
Fukuhara  (also  spelled  Hukuhara).  For  the  past  year,  Watanabe  has  implemented 
a  package  that  transforms  a  large  class  of  equations  into  P-functions.  He  has 
used  Kamke's  table  as  a  test  bed.  Several  examples  from  Kamke  using  this 
package,  are  given  in  Appendix  I. 


APPENDIX  I 


(D2 )  [DSK,  SWATAN] 

(C3)  showtime : trueS 
Time*  5  msec. 


( C 4 )  loadf i 1 e ( pmai n , f as! ) ; 

PMAIN  FASL  DSK  SWATAN  being  loaded 
Loading  done 
Time*  411  msec. 

(D4)  DONE 

{ C 5 )  batch( exampl , 2  ) ; 


Kanke  exanole  3?8a 


( C  6 ) 

Time 


(D6) 


k378a:x*(x-l)*(x+l)'2*,diff(y.x.2)+2*x*(x+l)*(x-3)*’diff(y.x)-2*( 
57  msec. 


2 

2  d  Y  dY 

(X  -  1)  X  (X  ♦  1)  —  ♦  2  (X  -  3)  X  (X  +  1)  -  2(X-1)Y* 

2  dX 

dX 


( C 7 )  lode2(k378a); 

dY 

2  (2  X  -  6)  -- 

d  Y  dX  2  Y 

we  solve  —  + - - - 

2  2  3  2 

dX  X-l  X  +  2  X  ♦  X 

SOLVE  FASL  DSK  MACSYM  being  loaded 
Loading  done 


0 


PREMAP  FASL  DSK  SWATAN  being  loaded 
Loading  done 

U 

we  use  Y  *  - 

2 

(X  ♦  1) 

dU 

2  --  2 
2  U  dX  d  U 

the  result  is  - -  - * - *  0 

2  X-l  2 

X  -  X  dX 


dY 

2  2  -- 

d  Y  dX  2  Y 

we  solve  —  - - ♦ - -  0 

2  X-l  2 

dx  x  -  x 

PEQPl  FASL  DSK  SWATAN  being  loaded 
Loading  done 


iP  POr.-TE  o,vv  •  r 


PHYPGM  FASL  DSK  SWATAN  being  loaded 
Loading  done 


x- 1 ) *ys0 


0 


(02)  [DSK .  SWATAN] 

(C3)  showtime : true$ 

Time*  4  msec. 

(C4)  batch( exampl , 7) ; 

( C5 )  K344:X^4*’DIFF(Y,X.2)+(EXP(2/X}-V*2)*Y=0; 

Time*  65  msec. 

2 

4  d  Y  2/X  2 

(05)  X  ---  +  (XE  -  V  )  Y 

2 

dX 


(C6)  loadfi1e(pmain,fasl ) ; 

PMAIN  FASL  DSK  SWATAN  being  loaded 


Loading  done 
Time*  413  msec. 

(06)  DONE 

( C 7 )  lode2( k344) ; 

2  2/X  2 

d  Y  (XE  -  V  )  Y 

we  solve  —  + - *  0 

2  4 

dX  X 


PIVTR  FASL  DSK  SWATAN  being  loaded 
Loading  done 

PTRINV  FASL  DSK  SWATAN  being  loaded 
Loading  done 
1 

we  use  T  «  - 
X 

dY 

2  2  ~ 

d  Y  dT  2  T  2 

the  result  Is  —  +  . - +  (XE  -  V  )  Y  *  0 

2  T 
dT 

dY 

2  2-- 

d  Y  dX  2  X  2 

we  solve  —  + - +  (XE  -  V  )  Y  •  0 

2  X 
dX 

PIVTR2  FASL  DSK  SWATAN  being  loaded 
Loading  done 

PTREXP  FASL  DSK  SWATAN  being  loaded 
Loading  done 


PDVTR  FASL  DSK  SWATAN  being  loaded 


the  type  is  hypergecmetri c 

the  solution  may  be  written  by  Riemann’s  P-functions  as  follows 
[0  1  INF  ] 

[  ] 

y=P  [  0  0  -  2  ]  (x) 

C  3 

[13-1] 

PNOPM  FASl  DSK  SWATAN  being  loaded 
Loading  done 


we  solve 

[0  1  INF  ] 

[  3 

y=P  C  i  o  -  i  3  (x) 

C  3 

[03-2J 
/  2 

[  (X  -  1) 

(K2  I  . . dX  +  Kl)  X 

3  2 

/  X 


continue?  type  y  or  n 

y; 

SIN  FASL  DSK  MACSYM  being  loaded 
Loading  done 

SININT  FASL  DSK  MACSYM  being  loaded 
Loading  done 

SCHATC  FASL  DSK  MACSYM  being  loaded 
Loading  done 
Time1  21253  msec. 

2 

(D7)  -  2  K2  X  LOG ( X )  +  K2  X  ♦  Kl  X  -  K2 

Time*  23672  msec. 

(D8)  BATCH  DONE 


(C9)  closefile(buffer , save) ; 


Loading  done 


SOLVE  FASL  DSK  MACSYM  being  loaded 
Loading  done 
U 

we  use  Y  *  - 
X 

2 

2  X  2  d  U 

the  result  is  U  (XE  -  V  )  ♦  —  *  0 

2 

dX 

2 

d  Y  2  X  2 

we  solve  —  +  (%E  -  V  )  Y  -  0 

2 

dX 

x 

we  use  T  *  XE 

dY 

2  --  2  2 
dY  dT  {V  -  T  )  Y 

the  result  is  —  +  --  -  -  *  0 

2  T  2 

dT  T 

dY 

2  --  2  2 
d  Y  dX  (X  -  V  )  Y 

we  solve  —  ♦  --  + - *  0 

2  X  2 

dX  X 

PCONFL  FASL  DSK  SWATAN  being  loaded 
Loading  done 

the  equation  is  confluent  type 

PCNEQO  FASL  DSK  SWATAN  being  loaded 
Loading  done 

the  solution  may  be  representable  by  Fukuhara's  P-function. 
[  INF  INF  0  ] 

[  ] 

[  1  ] 

[  -  XI  -  -  V  ] 

y-P  [  2  ]  (x) 

C  ] 

[  1  3 

[XI  -  V  ] 

[  2  ] 

PCFPTM  FASL  DSK  SWATAN  being  loaded 
Loading  done 
y»  Y  (X) 

B,  ABS ( V ) 

Time*  15604  msec. 

(02)  Y  (X) 

•  B.  ABS(V) 


[DSK.  SWATAN] 


{ C 3  )  showtime : trueS 
Time*  5  msec. 

{ C 4  )  batch; exampl , 12 ) ; 

( C 5  )  k406: 16* ( X*3-l }~2* - D I FF ( Y . X . 2 ) -2 ? * X • Y  =  0 ; 

Time*  41  msec. 

2 

3  2  d  Y 

(D5)  16  (X  -  1)  —  +  27  X  Y  *  0 

2 

dX 


Kamke  examrle 


( C6)  k406t:48«x*2»(x-l)~2*'diff(y,x12)+32*x*(x-l)~2*’diff(y,x)+9*x*y=0 
Time*  50  msec. 

2 

2  2  d  Y  2  dY 

(D6)  48  (X  -  1)  X  —  +  32  (X  -  1)  X--*9XY=0 


(C7)  loadfile(pmain.,fasl ) ; 

PMAIN  FASL  DSK  SWATAN  being  loaded 
Loading  done 
Time*  426  msec. 

( D7 )  DONE 

( C8 )  Iode2(k406t); 


w 

2 

d  Y 

we  solve  —  + 

dY 

2  -- 
dX 

- *f 

3  Y 

2 

3  X 

3 

i 

h 

dX 

16  X 

-  32  X 

& 

Cat 

SOLVE  FASL  DSK 

MACSYM 

bei  ng 

loaded 

Loading  done 

PHYPGM  FASL  DSK  SWATAN  being  loaded 
Loading  done 

the  type  is  hypergeometric 

the  solution  may  be  written  by  Riemann’s  P-functions  as  follows 
[0  1  INF  ] 

C  3 

[13  1  ] 

[-*--] 
y-P  [  3  4  3  ]  <  x ) 

C  ] 

[  1  3 

[  0  -  0  ] 

[  «  3 


PHGHP  FASL  DSK  SWATAN  being  leaded 
Loading  done 


PALGS  rASL  CSK 

S  WAT  AN 

being 

loaded 

Loading  done 
we  solve 

[  o  i 

INF 

] 

[ 

] 

[  1  1 

1 

] 

1/4 

[  -  * 

-  -- 

] 

y  *  (x  -  i)  p 

[  3  2 

12 

]  (x) 

[ 

] 

[ 

1 

] 

[  0  0 

- 

] 

[ 

4 

] 

PALG4  FASL  DSK 
Loading  done 

SWATAN 

being 

loaded 

poh-s:  .'asl  rsK 

Loading  done 

SWATAN 

being 

1 02ded 

1/4 

2 

K 1  T  (2  SQRT ( T 

♦  T  ♦ 

1)  +  SQRT ( 3 )  (T 

jr 

3  1/12 

( SQRT ( 3 )  (T  ♦  1)  -  2  SQRT ( T  ♦  T  ♦  1)) 


where  t  =  x~( 1/3) 

Time*  13449  msec. 

1/4  2 

K 1  T  (2  SQRT ( T  +  T  +  1)  ♦  SQRT( 3 )  (T 


K2  7  ( SQRT ( 3  )  (T  ♦  1)  -  2  SQRT ( T  ♦  T  +  1 ) ) 

3  1/12 

(T  -  1) 


Time*  15278  msec. 
(09) 


BATCH  DONE 


(CIO)  closef i le( buf f er , save ) ; 


(D2) 


[DSK,  SWATAN] 


( C 3  )  showtime : t“ue$ 
Time*  5  msec. 

( C 4 )  batch( exampl , 13) ; 


Kanke  exacole  ISO 


( C 5 )  K180 : X'2  * 
Time*  47  msec. 

(D5) 


DI F F ( Y,X .2)+3*X*  * D1FF{ Y,X. 1)  +  (X‘2  +  1-V‘2)*Y*F(X): 
2 

2  d  Y  dY  2  2 

X  ---  +  3  X  --  +  (X  -  V  +  1)  Y  =  F ( X ) 

2  dX 

dX 


( C 6 )  loadf i le( pmai n , fasl ) ; 

PMAIN  FASL  DSK  SWATAN  being  loaded 
Loading  done 
Time*  413  msec. 

(D6)  DONE 

( C 7 )  lode2( k 180) ; 

dY 

2  3--  2  2 

d  Y  dX  (X  -  V  +  1)  Y  F ( X ) 


we  solve  —  + - + -  - - 

2  X  2  2 

dX  XX 

dY 

23--  22 

d  Y  dX  (X  -  V  +  1)  Y 

first  we  solve  —  + - + -  =  0 

2  X  2 

dX  X 


SOLVE  FASL  DSK  MACSYM  being  loaded 
Loading  done 


PCONFL  FASL  DSK  SWATAN  being  loaded 
Loading  done 

the  equation  is  confluent  type 


PCNEQO  FASL  DSK  SWATAN  being  loaded 
Loading  done 

the  solution  may  be  representable  by  Fukuhara's  P-function. 

[  INF  INF  0  ] 

C  3 

[  3  ] 

[  -  XI  -  -  V  -  1  ] 

y-P  [  2  ]  (x) 

C  3 

[  3  ] 

[XI  -  V  -  1  ] 

[  2  ] 


PCFPTM  FASL  DSK  SWATAN  being  loaded 


Loadmq  done 

Y  (X) 

B.  ABS(V) 

y*  . 

x 


the  solution  of  the  homog.  eq.  is 


Y  (X) 

B.  ABS(V) 


PNONH  FASL  DSK  SWATAN  being  loaded 
Loading  done 

a  special  sol.  for  the  nonhom.  eq.  is 
X 

/ 

[ 

I  (%PI  F ( T )  J  (T)  Y  (X)  -  %PI  F ( T )  Y  (T)  J  (X))  dT 

]  ABS( V}  A8S { V)  ABS(V)  ABS(V) 

/ 

0 


continue?  type  y  or  n 

n; 

Time*  9943  msec. 

X 

/ 

[ 

I  (%PI  F(T)  J  (T)  Y  (X)  -  %PI  F(T)  Y  (T)  J 


(X))  dT 


Time*  12131  msec. 

(08) 

BATCH  DONE 

B.  ABS(V) 
+ - 

X 

(C9)  cl osef i 1 e ( buf f er , save  ) ; 

